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Inhomogeneous broadening of energy levels is one of the principal limiting factors for achieving 
"slow" or "stationary" light in solid state media by means of electromagnetically induced trans- 
parency (EIT), a quantum version of stimulated Raman adiabatic passage (STIRAP). Stark-shift- 
chirped rapid-adiabatic-passage (SCRAP) has been shown to be far less sensitive to inhomogeneous 
broadening than STIRAP, a population transfer technique to which it is closely related. We fur- 
ther optimise the pulses used in SCRAP to be even less sensitive to inhomogeneous broadening in 
a A-type three-level system. The optimised pulses perform at a higher fidelity than the standard 
gaussian pulses for a wide range of detunings (i.e. large inhomogeneous broadening). 

PACS numbers: 32.80.Qk, 33.20.Bx, 33.80.Be, 42.50.Hz, 42.50.Gy 



I. INTRODUCTION 

Atomic vapours have been extensively investigated for 
use as quantum information storage media by "slowing 
down" , or even "stopping" , a light pulse carrying quan- 
tum information. This is achieved by mapping the quan- 
tum state of the light to a long-lived spin state in an en- 
semble of atoms by means of a reversible process: electro- 
magnetically induced transparency (EIT) P, [2, H, 0, [1] , 
which is closely related to an adiabatic passage technique 
known as Stimulated Raman Adiabatic Passage (STI- 
RAP). The homogeneity of the atomic vapour atoms 
means that the light pulses used in EIT experiments 
in these systems can be tuned to specific transitions in 
the atoms and that these light fields will then interact 
strongly with all the atoms in the ensemble. 

Solid-state implementations of stored light by means 
of EIT may have a number of potential advantages over 
atomic vapour implementations. They have potential 
to greatly reduce (possibly eliminate) limitations on the 
storage lifetime, which is mainly due to atomic diffu- 
sion and Doppler velocities in atomic vapour systems. A 
higher atomic density has the potential to yield a stronger 
interaction between the light and atomic ensemble, whilst 
being more compact and simple to use/manufacture may 
indicate great scalability on the part of solid-state 
implementations. EIT and slow light have been demon- 
strated in rare-earth doped semiconductors e.g. Pr 
doped Y^SiO^ 0, i, whereas only EIT has thus far 
been shown in nitrogen-vacancy colour centers in dia- 
mond p. [loj. Solid-state media are, however, not with- 
out their own problems, the biggest of which is inhomo- 
geneous broadening of the energy levels, which leads to 
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a reduction in the number of atoms/centers that will be 
resonant with the incident optical field. This broadening 
occurs because each atom/center experiences a different 
electronic environment due to a number of factors e.g. 
the local strains in the crystal lattice. 

Recently a variation of STIRAP (the basis of 
EIT), namely Stark-shift-chirped rapid adiabatic passage 
(SCRAP), has been shown to perform coherent popula- 
tion transfer with a high fidelity for a ran ge of different 
detunings in a A-type three-level system That is, 

SCRAP is far more robust when there is large inhomoge- 
neous broadenings present in the ensemble. In SCRAP, a 
separate (pulsed) field is added to the STIRAP fields to 
induce Stark shifts in the energy levels and thus bring the 
required, initially off-resonant, transitions on resonance 
at specific times. We propose that a quantum version of 
SCRAP could surmount some of the limitations that in- 
homogeneous broadening places on "slow" light in solid 
state systems. We will examine this further in future 
work. 

In this work we make use of optimum control tech- 
niques similar to Khaneja et a. I (2005) [l^, in order 
to optimise the standard SCRAP pulses so as to min- 
imise the decrease in fidelity brought on by inhomoge- 
neous broadenings of the transitions. To measure this 
we simulate the SCRAP process population transfer be- 
tween two long lived ground states which experience large 
detunings due to inhomogeneous broadening. Our main 
result is that we can improve the average fidelity of pop- 
ulation transfer over a wide range of detunings for both 
the ground to excited state detuning and the ground to 
target state detuning (two-photon detuning). The opti- 
mal control pulses are thus tailored to provide effective 
state transfer (and thus EIT), in the presence of large 
inhomogeneous broadening. 

The paper is arranged as follows: In section |TT] the 
SCRAP technique is introduced. Section IIIII describes 
the optimisation methods used. Section IIVI shows the 
results for the optimised pulses. Conclusions are drawn 
in Section IVl 
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FIG. 1; The A-type three-level energy scheme. States |1) and 
1 2) are coupled by the pump pulse P which has a detuning A p 
from being exactly on resonance. Similarly states |2) and |3) 
are coupled by the Stokes pulse S which has a detuning As- 
State 1 2) is short lived with spontaneous emission occurring 
out of the system. 



II. SCRAP 

The Stark-chirped rapid adiabatic passage technique 
(SCRAP) was first proposed by Yatsenko et al. and 
implemented by Rickes et al. [Ol, in two- level systems 
as an efficient method for complete population transfer 
between two states. It was later shown by Rangelov et al. 
(2005) fu\, that a Stark shifting pulse can also be used to 
achieve complete population transfer through adiabatic 
passage in a three-level system, thus providing an alter- 
native to STIRAP. 

The A-type three-level system (Figure [1]), is comprised 
of two long lived ground states, one the initially pop- 
ulated state (|1)), and the other the target state (|3)), 
and an excited state (|2)). In general the excited state 
has a short life time: STIRAP is effective at complete 
population transfer since it avoids populating the excited 
state by having the system evolve along a dark-state. In 
both SCRAP and STIRAP states |1) and |2) are coupled 
by the "pump" laser pulse whilst states |3) and |2) are 
coupled by the "Stokes" laser pulse. The frequencies of 
these two classical laser fields are typically not exactly 
on resonance with their respective transitions, and have 
detunings Ap and A5 for the pump and Stokes lasers 
respectively. Herein lies the advantage that SCRAP has 
over STIRAP: STIRAP requires exact two-photon reso- 
nance (Ap — As = 0) in order to be effective, whereas 
SCRAP has a larger tolerance for two-photon detuning. 
In SCRAP a third strong far-off-resonance laser pulse, 
the Stark pulse, is also applied. The Stark pulse induces 
a Stark shift in the energy of the excited state (essen- 



tially leaving the energy of the lower levels unchanged), 
thus bringing the initially off-resonant transitions into 
resonance. The Stark shift causes the diabatic energy 
of state |2) to cross those of states |1) and |3), allowing 
population transfer from |1) to |2) and then to |3), com- 
pleting the population transfer. It is thus clear that the 
diabatic energies of state |1) and |2) must cross before 
the crossing of energies of states |2) and |3). It is also 
clear that decay out of the excited state will play a role 
in SCRAP, as opposed to STIRAP. 

In the rotating wave approximation, the Hamiltonian 
of the A-type three-state system depicted in Figure [T] is 
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where flp(t) and ^sit) £^re respectively the pump and 
Stokes laser field Rabi frequencies and 5*2 (t) is the Stark 
shift in the energy of the excited energy level |2) due to 
a third far-off-resonant laser pulse (Stark pulse). For the 
examples here it is assumed that the Stark shift is nega- 
tive, 5*2 (t) < 0. The detunings of the pump and Stokes 
fields are Ap and A5 respectively. The imaginary term 
iT describes the losses from |2) due to spontaneous ra- 
diative decay out of the three-level system. The Stark 
shifts of the energy levels for states |1) and |3) are not 
included here because Stark shifts in ground and meta- 
stable states tend to be much smaller than those of ex- 
cited states. 

In Rangelov et al. (2005) till]), gaussian pulse shapes 
were used for all the pulses, with identical peak values, 
ilo, for the Rabi frequencies of the pump and Stokes 
pulses, 
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The peak of the Stark pulse (maximum Stark shift of Sq) 
is taken to be at i = 0, and as such the pump and Stokes 
pulses peak at times Tp and Tg respectively. The pulse 
durations are determined by Tp, T5 and Tst, where it 
was taken that the pump and Stokes pulses have equal 
duration Tp = Ts and the Stark pulse has twice their 
duration Tst = 2Tp. The unit of time was defined as 
Tp and the unit of frequency as 1/Tp. These gaussian 
pulses served as exemplars for the initial pulses used in 
our optimising routine, see section [Till 
The diabatic energies of the states are 



EiD = 0, 
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which with the condition that the Stark shift is negative 
(52 (t) < 0), dictates that the diabatic energy of state 
2), -E|2), will only cross those of states |1) and |3) when 



and 



5o > Ap > 0, 



So > As > 0. 
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With the condition that the Stark shift is negative two 
distinct situations can arise: the two-photon detuning 
can be negative 



Ap > > (Ap - As) > Ap - 5o, (10) 



or positive 

Ap > (Ap - As) > > Ap - S-Q. (11) 

In the examples used here to explain the SCRAP tech- 
nique, and for the pulse optimisation, only the first case 
(llOp . where the two-photon detuning is negative will be 
presented (Ap < As). In the case that the two-photon 
detuning is positive (jlip , the order of the pulses in stan- 
dard SCRAP must be run in reverse to what will be 
shown here, 



Our goal is to achieve efhcient state transfer for as 
wide a variety of detunings as possible, for which SCRAP 
is the ideal technique. The population transfer for the 
ideal situation with no detuning, no decay, and using 
pulse parameters out of [ll| is shown in Figure [2l In 
such a situation with no detuning STIRAP would be the 
preferred technique. For a more detailed explanation of 
three state SCRAP please refer to section III in Rangelov 
et al. (2005) 



III. OPTIMISATION THROUGH OPTIMAL 
CONTROL 

The above mentioned SCRAP process has the advan- 
tage that it can tolerate large two-photon detunings, but 
unlike STIRAP it suffers from decay from state |2). In 
this section we will use an optimisation technique, based 
on the GRAPE algorithm by Khaneja et al. (2005) [H, 
to optimise the transfer fidelity for a set of detunings. 
The state of the three level system is characterised by 
the density operator p (t) with the Liouville-von Neuman 
equation of motion 
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Where Hq is the free evolution Hamiltonian and the 
Hk are the Hamiltonians corresponding to the m con- 
trol fields (the Stokes, pump and Stark fields). u{t) = 
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FIG. 2: The evolution of the population in the 3-state-system 
(state |1) solid line, |2) dashed line and |3) dash-dotted line) 
where the pulses used were the original gaussians with the fol- 
lowing parameters: Ap = 30/Tp, As = 45/rp, So = 200/rp, 
flo = 50/rp, Ts = Tp, Tst = 2Tp, = -Tp, = -2rp. 
There was no decay out of state |2), (F = 0). 



(mi [t) ,U2{t) , . . . , Um [t)) is the vector of control ampli- 
tudes. We discretised the transfer time T into N steps 
of length Ai = T/N , and assume that the amplitude 
for each control field is constant during each time step. 
Instead of optimising these amplitudes directly, as de- 
scribed in |12| . we define each amplitude vector in terms 
of q gaussians. 



(13) 



that sum to create the pulse for the specific control field, 
and optimise the parameters of these gaussians. The aim 
is to find the parameters (ft.„,fe, T„_fc, an,k) that will, given 
the initial density operator p (0) = po, maximise the over- 
lap of the density operator after a time T, p(T), with a 
target density operator C . The overlap is measured by 
the standard inner product, thus the performance index 
$0 is given by 



^n = {C\p{t)). 



(14) 



During each time step j the evolution of the system is 
given by the propagator 



Uj — exp 



iAt [Ha + ^Uk (j) Hk 
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(15) 



The performance index can then be written as 
*o - (C/J+i . . . uIcUn ■ ■ ■ C/,+1 \U,... UipoUl . . . U]). 
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From [1^ we have that 
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but we are interested in the gradient with regard to hn^k, 
T„ i- and an,k- 
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and similarly 
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The performance $0 increases if we choose 



hn,k = hn.k + £ 



(21) 



with e a small step size, and similarly for r„ ^ and (T„jt. 

We made use of a Matlab routine 
"minFuncBC", written by Mark Schmidt 
(http : //www. cs . ubc . ca/~schmidtm/Sof tware/minFunc 
to perform the final optimisation step pip. It makes 
use of a quasi-Newton method, the BFGS method, 
and accepted as input $0 and the vector of deriva- 
tives with respect to the optimising parameters, 

(5<E>o li^o iS^o 
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TABLE I: The parameters for each of the 9 gaussians that 
constitute each of the pulses. The amplitudes of the gaussians 
for a given pulse were all eaqual: h„^p = h„^s = 0.23r2o and 

hn^St ~ O.23S0. 
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FIG. 3: The original gaussian pulses (solid lines) with So = 
200/Tp, Qo = 50/Tp, Ts = Tp, Tst = 2Tp, Tp = -Tp, 
Ts = —2Tp. The dashed lines are the pulses optimised for the 
detunings shown in Figure [5l 



no additional points result in better pulses. The total 
h'Jiaiformance function used for evaluating the efficiency 
of the pulses during the optimisation on a number d of 
detunings was taken as the average of the performance 
functions (fT4| for each of the d detuning points. 
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IV. RESULTS 

In order to optimise the pulses for as large a detuning 
space as possible, an initial single point in the detun- 
ing space was chosen by searching for the detuning point 
with optimised pulses that performed the best over the 
whole chosen detuning space. This entailed calculating 
the fidelity, using equation , for each point in the de- 
tuning space we would like to optimise, and then taking 
the average. Once the first point was found a second 
point would be chosen by searching for the second point 
that, together with the first, would result in the best 
optimised pulses. This process would be repeated until 



For the initial pulses used in our optimal control rou- 
tine we used 9 gaussians of equal amplitude to approxi- 
mate each of the original pulse shapes ^ and Q. 
That is, in equation (fT3|) . q = 9 and k S {P,S,St} 
(the probe, Stokes and Stark fields). The parameters 
of these gaussians are given in Table HI The original 
pulses (solid lines in Figure [3]) had the following param- 
eters: So = 200/rp, Qo = 50/Tp, Ts = Tp, Tst = 2Tp, 
Tp = -Tp, Ts = -2Tp, Tst = 0. 

The maximum of the optimised pulses was constrained 
to that of the original pulses {So — 200/Tp, flo — 
50/Tp). If this is not done pulses can be made arbi- 
trarily efficient by simply increasing the maximum Rabi 
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frequency. This is not physically possible due to adia- 
baticity constraints. Furthermore the optimisation was 
constrained to prevent the optimised pulses from becom- 
ing too narrow by setting the lower bound on the width 
of the pulses to 50% of the original pulse width. The 
optimised pulses are shown as the dashed lines in Figure 

m 

The optimised pulses were then used to evaluate the 
fidelity of the population transfer from state |1) to state 
3) for a range of detunings, keeping in mind that the 
chosen pulse ordering restricts the choice of detunings 
(fTU)l . As can be seen from Figures|3]and[5]the total "area" 
of detunings where population transfer is at all possible 
(the parts that are not black) is greatly increased when 
using the optimised pulses. The efficiency of the pulses 
can be gauged in a number of ways: 

1. measuring the normalised "area" in detuning space 
where the pulses result in a fidelity greater than 0.8 
for population transfer. For the original SCRAP 
pulses A'^Q g = 0.131, whilst for the optimised 
pulses = 0.178 was obtained, an increase of 
35.5%. 

2. measuring the average fidelity over the whole de- 
tuning space, i.e. the sum of the fidelities for each 
point, divided by the number of points. For the 
original SCRAP pulses F°^' = 0.219, whilst for 
the optimised pulses F°p* = 0.321, an increase of 
46.6%. 

3. in Figure [6] a plot of the percentage increase for 
each point in the detuning space is presented. The 
logj^Q of the percentage increase is used since some 
points had an initial fidelity of virtually zero and 
ended with a percentage increase of « 10^^%. 
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FIG. 4: The fidelity for a range of detunings, using the original 
gaussian pulses for SCRAP (shown as the solid lines in Figure 
O, with So = 200/Tp, fto = 50/rp, Ts =Tp, Tst = 2Tp, 
Tp = — Tp, Ts = —2Tp. The decay out of state |2) was F = 
l/Tp. 
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FIG. 5: The fidelity for a range of detunings, using SCRAP 
pulses optimised for the detunings indicated (white cirlcles). 
Optimised pulses shown as the dashed lines in Figure |3] The 
decay out of state |2) was F = 1/Tp. 
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FIG. 6: The logj^Q of the percentage increase in fidelity be- 
tween standard SCRAP (Figure \Q and optimised SCRAP 
(Figure [5]) for each point in the detuning space. 



V. CONCLUSIONS 



It is clear that Stark-shift Chirped Rapid Adiabatic 
Passage (SCRAP) is a useful process to employ when 
trying to overcome inhomogeneous broadening of energy 
levels in a system undergoing state transfer. We have 
shown that the standard SCRAP pulses can be optimised 
so that a larger inhomogeneous broadening can be com- 
pensated for and so that the overall fidelity of state trans- 
fer for a range of detunings will be increased. Future work 
will aim to translate SCRAP into the quantum domain 
where the pump pulse is replaced by a quantum probe 
field carrying quantum information that is to be stored 
in the spin coherence of the atomic system. 
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